Environment-independent decoherence rate in 
classically chaotic systems. 

Rodolfo A. Jalabert 1 and Horacio M. Pastawski 1 ' 2 
l Institut de Physique et Chimie des Materiaux de Strasbourg, 
UMR 7504, CNRS-ULP, 23 rue du Loess, 67037 Strasbourg Cedex, France 
2 Facultad de Matemdtica, Astronomia y Fisica, Universidad 
Nacional de Cordoba, Ciudad Universitaria, 5000 Cordoba, Argentina 

February 1, 2008 

Abstract 

We study the decoherence of a one-particle system, whose classical 
correpondent is chaotic, when it evolves coupled to a weak quenched en- 
vironment. This is done by analytical evaluation of the Loschmidt Echo, 
(i.e. the revival of a localized density excitation upon reversal of its time 
evolution), in presence of the perturbation. We predict an exponential 
decay for the Loschmidt Echo with a (decoherence) rate which is asymp- 
totically given by the mean Lyapunov exponent of the classical system, 
and therefore independent of the perturbation strength, within a given 
range of strengths. Our results are consistent with recent experiments of 
Polarization Echoes in nuclear magnetic resonance and preliminary nu- 
merical simulations. PACs: 03.65.Sq, 05.45. +b, 05.45. Mt, 03.67.-a 

The coupling of a system to environmental degrees of freedom plays an im- 
portant role in many areas of physics. Already on a classical level, it leads to 
fluctuations, damping and irreversibility. In quantum mechanics the environ- 
mental coupling induces decoherence, destroying quantum superpositions and 
reducing pure states to a mixture of states Jj[] . It is then not surprising that the 
concepts of environment, decoherence and irreversibility have been the object 
of scholar discussions for long time Renewed interest has been fostered by 
the crucial role that decoherence plays in the problem of quantum computation 
[PJ and by the technical advances that make it possible to perform experiments 
envisioned as gedanken for long time. 

Experiments with Rydberg atoms in a microwave cavity Q] allow to observe 
the progressive decoherence in a quantum measurement problem, while analysis 
of conductances through semiconductor microstructures Q make it possible to 
address the "which path" problem in a solid state environment. In addition, 
nuclear magnetic resonance (NMR) offers unlimited possibilities for the study 
of decoherence and irreversibility in a tailored environment. The phenomenon 
of spin echo, shows how an individual spin, in an ensemble, loses its "phase 
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memory" as a consequence of the interaction with other spins that act as 
an "environment" . The failure of recovering the initial ordered state in a time 
scale Tii manifest the effect of a complex (many-spin) environment on the re- 
versibility of simple systems. A further conceptual breakthrough was enabled 
by experiments that revert JjJ and control^ the whole entangled state of the 
strongly interacting nuclear spins to obtain the NMR Polarization Echo. We 
will briefly discuss the physics involved. 

A local spin excitation \tp) created at time t — spreads out through the 
crystal under the action of a many-spin Hamiltonian Ho allowing exchange 
between spins. This complex quantum evolution is macroscopically assimilated 
to a "spin-diffusion" Q process (consistently with the usual hypothesis of micro- 
scopic chaos describing many particle systems ]l0|). At time i,a radio- frequency 
pulse sequence produces a new effective Hamiltonian — (Ho + £). Here S, a 
perturbation containing the pulse imperfections and residual interactions with 
additional spins, can be made very smallQ. Hence, after the pulse at t, one gets 
an implementation of the gedanken backwards dynamics proposed by Loschmidt 
in his argument against the Boltzmann's H-theorem. At time 2t, one measures 
a maximum in the return probability to the initial state that we call Loschmidt 
Echo (LE), 

M(t) = (tJj\ e ^o+v)t/h e -m t/h ^ 2 ^ 

The build up of the LE depends on a precise interference between the "diffusive" 
wave-packets e - m ° t/h \ip) and e -^+^)t/h ^ ^ that ig degraded by S. Clearly, 
M(t) should be a decreasing function of the elapsed time t before the reversal 
of Ho with a decoherence rate < I/T2. A surprising outcome of the 

experiment Q is that, for small S's, l/r^ only depends on the intrinsic properties 
of the system (that is, on Ho)- 

In this work, we develop a simple analytical model exhibiting the indepen- 
dence of the decoherence rate on the perturbation found in the experiment. 
The system is represented by a single-particle Hamiltonian Hq whose underly- 
ing classical dynamics is strongly chaotic. This is clearly an oversimplification 
respect to the many-body Hamiltonian of interest, but still it introduces enough 
complexity in the intrinsic system (quantified by its mean Lyapunov exponent 
A) which is absent in simpler dissipative systems, where Ho is integrable. Plac- 
ing ourselves between the limits of a trivial and a many-body Ho allows us to 
have a tractable model and explore the influence of classical chaos in quantum 
dynamics. To account for "non-inverted" part of the Hamiltonian evolution we 
consider an Hermitian operator E representing the coupling with a quenched 
environment acting in the backward evolution (from t to 2t). This approach is 
not only consistent with the experimental situation but it is also able to provide 
a new insight into the problem of decoherence because the calculation can be 
handled within the precise framework of the Schrodinger equation. In contrast, 
most of the previous studies of decoherence use extremely simple Hamiltonian 
systems (l| interacting with a dissipative environment (e.g. stochastic noise) 
which justifies the use of a master equation for the reduced density matrix. In 
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this context, by discussing the entropy growth of a dissipative system, Zurek and 
Paz hinted at the importance of the chaotic classical dynamics in setting the 
characteristic time scales for decoherence. This is consistent, under conditions 
that we specify bellow, with our results for the time decay of M(t). However, 
since we use a purely Hamiltonian approach, our conceptual framework is very 
different. 

As in the experiment, we start with a localized state in a <i-dimensional 
space, 



= o) = (_§*) 



. d/4 



exp 



ip • (r - r ) - 5^ (r - r ) 



(2) 



centered at ro, with dispersion ex. The momentum po selects the energy range 
of the excitation. This choice also renders the calculations tractable. The time 
evolution of the initial state is best described using the propagator K(r, r; t) = 
(r| e - mt / h |r) by 



V(r;t) = / df]f(r,f;t) ^(r;0) . (3) 

Using the Hamiltonian Ho + E or TCq in the propagator K yields tpHo+'s or tp^ , 
respectively. We take E as a static disordered potential given by Ni impurities 
with a Gaussian potential characterized by the correlation length £, 



S = V (r) - L J^ym -P [-M^) 2 \ ■ ( 4 ) 

The independent impurities are uniformly distributed with density n; = A^/V, 
(V is the sample volume). The strengths u a obey {u a Uf}) = u 2 S a f3- This assump- 
tions simplifies analytical evaluation of the ensemble average of the observable 
M(t). We stress that we are not simply describing the physics of disordered 
systems (which is obviously phase coherent), since the potential V'(r) acts in 
the backwards propagation but not in the forward path. 

As a calculational tool, we use the semiclassical approximation for K(r, r; t), 
as the sum over all the classical trajectories s(r, r;t) joining the points r and r 
in a time t @, [if]; 

K(r,f;t) = J2 K s(r,f;t) , with (5) 

s(r,r;t) 

K s (v,f;t) = (^y^C^exp [^S s (v, r; t) - f p,] , 

valid in the limit of large energies for which the de Broglie wave-length (Xp = 
2n/k = 2Tth/po) is the minimal length scale. iS is the Hamilton principal func- 
tion (action) specified by the integral of the Lagrangian S s (r,f;t) — J Q diC 
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along the classical path, /i is the Maslov index that counts the number of con- 
jugate points (focal points), and the Jacobian C s = |det23 s | accounting for the 
conservation of the classical probability, is expressed in terms of the initial and 
final position components j and i as (B s )ij = —d 2 S s /dridfj. This approxi- 
mation gives the wave function with great accuracy up to very long times H|. 
Besides, it provides the leading-order corrections in h due to V(r), in the limit 
of fc£ 3> 1, from the classical perturbation theory for the actions [^6, 17 
Using the initial wave-function of Eq. (||) and Eq. (||) we can write 



^(r; t) = (47ra 2 ) d/4 £ K s (r, r ; t) exp (p, - Vo f 

»(r,r ;i) 



(6) 



where we used that dS / dr.i\- =ra — —pi, and neglected the second order terms 
of S in (r — r ). This is justified under the assumption that £ ^> a ^> Xp, i.e. 
an initial wave-packet concentrated in a smaller scale than the fluctuations of 
V(r). In Eq. (^), only trajectories with initial momentum p s closer than h/a 
to po are relevant for the propagation of the wave-packet. 
The semiclassical approximation to the LE is: 



exp 



2h< 



dr£cy 2 C, 1/2 exp 



((p s - Po) 2 + (ps - Po)' 



- (S s - S s ) - — (n s 



(7) 



and involves two spatial integrations and four trajectories. 

The perfect echo of E = is already obtained considering only trajectories 
s = s which leaves aside terms with a highly oscillating phase: 

2 



M E=0 (t) = 



dr^C s exp -g^(p s -p ) 



= 1. 



(8) 



The integration requires the change from final position variable r to initial 
momentum p using the Jacobian C. 

In the coupled case the square modulus requires a the second integration 
variable r'. We see that only the terms with slightly perturbed trajectories 
s = s (as well as s' = §') survive the average over impurities. Thus, 



M(t) ~ (^) d / dr / dr' ^C s C s ,exp 



exp 



(< 



s,s' 

Ps - Po) 2 + (Ps' - Po) 2 



-(AS S -AS S ,) 
n 



(9) 



where = — J Q dt V(q s (t)) and AS S > are the phase differences, along the 
trajectories s and s', resulting from the perturbation V, From Eq. (^|) we see 
that we can decompose M into 

M(t) = M ad (t) + M d (t) , (10) 
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where the first term (non-diagonal) contains trajectories s and s' exploring 
different regions of phase space, while in the second (diagonal) s' remains close 
to s. 

In the non- diagonal term the impurity average can be done independently 
for s and s' . For trajectories longer than £ the phase accumulation AS S results 
from uncorrelated contributions, and therefore can be assumed to be Gaussian 
distributed jl^, [l7| . The disorder contribution involved in Eq. (|^) is then given 

by 



(exp 



n 



exp 



1 

" 2h 2 



dt 



di' Cy(\q s (t) - q s (t')\) 



where the correlation function of the disordered potential is 



^(|q-q'|) = <y(q)nq')) = 



( 47r £2)d/2 



exp 



-^(q-q') S 



(11) 
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The change of variables q — vt and q' = vP yields two integrals along the 
trajectory s. Since the length L s of the trajectory is supposed to be much larger 
than £, the integral over q — q' can be taken from — oo to +oo, while the integral 
on (q+q')/2 gives a factor of L s . Assuming that the velocity along the trajectory 
remains almost unchanged respect to its initial value i>o = po/m = L s /t , one 
gets, 



M 



ml 



drY^Cs 



exp 



W (Pa - Po)' 



exp 
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exp 



-tv /i 



(13) 



In analogy with disordered systems |16|, |Tj we have defined the typical length 
over which the quantum phase is modified by the perturbation as: 



fc2 2 



dq C(q) 



u 2 rii 



(14) 



We then see that M ad (t) has its time scale determined by E (through I), 
In computing the diagonal term M d (t) we use the expansion 



AS S - AS S , = / df VV(q s {t)) ■ (q s (t) - qy (t)) 



(15) 



since the trajectories s and s' remain close to each other. The difference between 
the intermediate points of both trajectories can be expressed using B : 



q s (f) - qy(t) = B~ 1 (t)(p s - p s /) = B-^B^t) (r - r') . 



(16) 
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In the chaotic case the behavior of S _1 (f) is dominated by the largest eigenvalue 
e xt . Therefore we make the simplification B^ 1 it)B{t) = exp [X(t — t)] I, where 
/ is the unit matrix and A the mean Lyapunov exponent. Here, we use our 
hypothesis of strong chaos which excludes marginally stable regions jD| with 
anomalous time behavior. Assuming a Gaussian distribution for the extra phase 
Of Eq. (ff| 



(exp 



~(AS S -AS S ,) 



exp[- 27-7 f dt [ At' exp [X(t + t' - 2b)] 
Jo Jo 



cwCM*)-^)!) ( r - r 'r 



We are now led to consider the "force correlator" 
C vt ,(|q-q'|) = (W(cO-W(q')) 



(4^2^/2 



2e 



exp 



'\2 



(17) 



(18) 



We change from the variables t and V to the coordinates q and q' along the 
trajectory s and use the fact that C V y is short-ranged (in the scale of £) to 
write, 



M c 



(t) ~ (^) y dr y dr' J] C 2 exp [-^ (p s - p ) 2 ] exp [-^ (r 



(19) 



where ^4 = (d— l)M 2 n i /(4Awo(47r^ 2 )^ -1 ^ 2 ) results from the f and t' integrations 
of Eq. ([ll]) in the limit Xt ^> 1. The last line comes from Gaussian integration 
over (r — r'). The factor C 2 reduces to C s when we make the change of variables 
from r to p. In the long-time limit C7 1 cx e At , while for short times C.T 1 = t/m. 
Using a form that interpolates between these two limits we have 



Af d (i) dp (^f) d ' 2 exp [-At]} exp [-#(p - p ) 2 

~ ^4exp [—At] , 



(20) 



with A = m/{A d /H). Since the integral over p is concentrated around po, 
the exponent A is considered constant . The coupling S appears only in the 
prefactor (through A) and therefore its detailed description is not crucial in 
discussing the time dependence of M d . The t factor in A induces a divergence 
for small t. However, our calculations are only valid in the limit Xt 1. Long 
times, (of the order of the Ehrenfest time tg = A -1 ln[fca] where a is a length 
characterizing Tto), are also excluded from our analysis since we run into the 
failure of the diagonal approximation. 
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Our semiclassical approach made it possible to estimate the two contribu- 
tions of Eq. @ to M(t). The non-diagonal component M nd (t) is the dominant 
contribution in the limit of small S. In particular, it makes M^ = o(t) = 1 in 
Eq. (||). The small values of £ are not properly treated in the semiclassical 
calculation of the diagonal term M d (t). While increasing the coupling £ the 
crossover from M nd to M d is achieved when / becomes smaller than vq/X. This 
condition is compatible with the assumption that, in the limit fc£ ^> 1, clas- 
sical trajectories shorter than the perturbation's "transport mean-free-path" 
l tr . = A(k£) 2 l are not affected [|l6| by the quenched disorder. For strong E 
the perturbative treatment of the actions is also expected to break down. 

We can now establish our main conclusion. In a system that classically 
exhibits strong chaos and can be characterized by a mean Lyapunov exponent 
A, a small random static perturbation may destroy our control of the quantum 
phase at a rate, 



— = - lim - lnM(<) = A , (21) 

Tijj t^oo t 

provided that the time is taken in the interval A -1 <C t <C t^, the perturbation 
presents long-range potential fluctuations (fc£ ^ 1) and a strength quantically 
strong (I <C vq/X) but classically weak (v /X <C Z t i ) • Notice that the thermody- 
namic limit, V — > oo, is required to take t arbitrarily large. 

The various restrictions for the validity of our result provide stringent con- 
ditions for its numerical verification. A disordered system represented by a 
tight-binding model with the topology of a torus exhibits a characteristic time 
of decay of M(t) depending on the disorder (Ho), but not on £ (giving the 
change of the magnetic flux piercing the torus) |l8) . Though subject to finite 
size limitations, the results show an environment independent behavior when the 
perturbation exceeds a critical value defined by the parametric correlations of 
the spectra jl9) . Even if our calculations and these preliminary numerical stud- 
ies deal with single-particle Hamiltonian, we expect that the generic behavior 
that we found is robust when considering Hamiltonians with larger complexity, 
like the many-particle case which is most relevant to the experiments motivating 
our work. 

The field of Quantum Chaos deals with signatures of the classical chaos on 
quantum properties. The most widely studied properties have been the spectral 
correlations |2(if , the wave function scars and the parametric correlations 
[ fl9| . In contrast, the studies of the temporal domain have been less developed, 
mainly because of the lack of clear quantities as those of steady state |22j. The 
work of Peres [^3| was a partial success in that direction. He distinguished 
regular and irregular dynamics on the basis of the asymptotic properties of a 
perturbation dependent overlap when V is finite and applied it to simple sys- 
tems. While there have been further attempts to define a dynamical sensitivity 
to perturbations in other systems and observables p4| ] , we are not aware of other 
predictions of a manifestation of the classical Lyapunov exponent in a Hamilto- 
nian system, as we did with Eq. (plj). On view of this result, we think that the 
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issue of dccohcrcnce by quantum evolution in classically chaotic systems, both 
with strong chaos and with marginal stability, deserves a more thorough exam- 
ination. Studies with other analytical and numerical techniques should clarify, 
among other aspects, the effects of different specific perturbations, the subtle 
effects of thermodynamic limits, the corrections due to Anderson localization, 
and the different temporal laws observed in one-body and many-body systems. 
This understanding of dynamical manifestations of chaos in the quantum world 
is decisive in the efforts to limit the experimental effects of decoherence and 
irreversibility. 

The authors have benefited from many fruitful discussions with O. Bohi- 
gas, P. Leboeuf, C. Lewenkopf, J. P. Paz and K. Richter. HMP is affiliated 
with CON1CET. We also acknowledge the financial support of the French- 
Argentinian program ECOS-Sud. 
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